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I.  BACKGROUND 


The  work  detailed  in  this  report  is  part  of  an  ongoing  effort  to  develop  models 
predicting  the  incidence  of  decompression  sickness  (DCS)  (1*5).  Report  I  of  the  series 
postulated  DCS  as  being  a  random  event  subject  to  the  laws  of  probability,  in  contrast  to 
the  traditional  view  of  DCS  as  a  deterministic  event  Semi-theoretical  models 
(describing  the  exchange  of  inert  gas  between  blood  and  tissue)  were  introduced  to 
predict  the  probability  of  decompression  sickness  [P(DCS)]  for  any  given  dive  profile, 
with  the  histories  of  depth  and  breathing  gas  composition  considered  the  sole 
independent  variables.  A  %rther  important  innovation  of  Report  I  was  the  use  of 
cumulative  risk  integrals  as  determinants  of  total  P(DCS)  in  the  models.  In  other  words, 
P(DCS)  was  considered  to  increase  over  the  course  of  the  dive  as  risk  accumulated 
according  to  the  model's  "rules".  These  models  were  fitted  to  over  1700 
well-documented  air  dives. 

In  Report  n  the  most  successful  of  the  models  from  Report  I  was  used  as  a 
predictor  of  P(DCS)  to  generate  recommended  air  diving  tables.  In  Report  In,  this 
same  model  was  used  to  compare  the  P(DCS)  of  dives  In  the  current  U.S.  Navy,  British 
Royal  Navy,  and  Canadian  Forces  air  tables.  In  Report  IV,  the  data-fitting  and  analysis 
were  extended  to  "saturation"  diveS,  and  it  was  demonstrated  that  the  modelling 
approach  outlined  in  Report  I  could  satisfactorily  predict  the  risks  of  dives  ranging  from 
less  than  a  minute  to  more  than  a  day  in  duration.  Accordingly,  recommended  tables  for 
saturation  dives  were  included  with  this  report.  Report  V  is  an  examination  of  an 


alternative  set  of  risk  models,  in  which  risk  is  viewed  as  something  that  depends  on  a 
single  incident  of  excess  dissolved  gas  at  some  instant  during  the  dive.  These  models 
were  found  to  be  less  adept  than  ours  at  fitting  heterogeneous  data  sets  such  as  those 
assembled  for  Report  IV. 

In  this  report  we  compare  repeat  and  multilevel  dives  with  single  dives,  asking 
whether  a  single  model  can  adequately  describe  the  various  sorts  of  dives.  A  repeat  dive 
is  defined  here  as  a  series  of  two  or  more  descents  separated  by  an  interval  of  less  than 
12  h  on  the  surface.  We  define  a  multilevel  dive  as  a  series  of  two  or  more  descents 
separated  by  an  interval  of  less  than  12  h  at  a  shallow  depth.  All  of  the  data  assembled 
for  this  report  consist  of  wet  dives  on  oxygen/nitrogen  performed  in  U.S.  Navy,  Canadian 
Forces,  or  British  Royal  Navy  trials.  Various  nitrogen/oxygen  breathing  mixes  were 
used,  although  most  were  air  dives  and  most  of  the  remainder  were  on  0.7  atm  oxygen, 
the  balance  nitrogen  (henceforth  to  be  called  simply  ”0.7").  Only  data  reported  since 
1978  were  used,  in  accordance  with  our  perception  that  possible  DCS  symptoms  are 
judged  by  substantially  different  standards  now  than  in  the  past.  Therefore,  we  did  not 
use  the  data  from  reports  I  or  IV  of  this  series  because  they  all  predated  1970. 

Since  this  study  includes  dives  on  several  gas  mixes,  it  became  necessary  to  test  the 
validity  of  our  conception  of  how  P(DCS)  is  affected  by  the  composition  of  the  breathing 
gas  mixture.  This  was  done  by  comparing  dives  on  different  breathing  gas  mixes  and 
asking  whether  data  collected  using  various  gas  mixes  can  be  described  by  a  single 
predictive  model. 
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The  comparisons  mentioned  above  were  perlonned  u  ng  objective  statistical  tests 
after  each  available  data  set,  or  combination  of  data  sets,  had  been  described  by  fitting 
the  probabilistic  models  to  it  The  indicated  results  are  that  the  models  do  not 
discriminate  between  single  and  repeat  dives  or  between  dives  on  different 
nitrogen/Ciorgen  breathing  mixes,  t'"  t  is,  it  does  not  ~>^em  necessary  to  use  different 
values  of  the  fitted  parameters  when  moving  from  one  category  of  dives  to  another. 
However,  none  of  the  models  used  in  this  report  can  adequately  describe  our  only  set  of 
multilevel  dives  while  simultaneously  providing  an  adequate  description  of  any  of  the 
other  data  s<ets.  Therefore,  we  carmot  say  with  confidence  that  the  present  models  are 
suitable  for  predicting  the  outcomes  of  multilevel  dives. 
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IL  MATHEMATICAL  MODELS 


The  models  have  been  described  in  previous  reports;  only  a  short  review  will  be 
presented  here.  In  the  simplest  conceivable  model,  every  dive  in  a  data  set  has  the  same 
P(DCS),  and  the  reasonable  course  is  to  set  this  value  of  P(DCS)  equal  to  the  observed 
average  P(DCS)  for  the  data  set  We  call  this  one*prj^meter  model  the  ”nuil  model" 
and  use  it  as  a  minimum  standard  that  any  creditable  semi-theoretical  model  must 
out-perform. 

The  models  of  real  interest  are  those  that  include  descriptions  of  gas  exchange 
kinetics.  In  such  a  model,  evaluation  of  the  safety  of  a  dive  is  accomplished  b>'  relating 
the  entire  dive  prefUe  to  the  probability  of  DCS  by  a  "risk  function": 

DCS  «  1.0  -  exp(  -  /q*  r  dt )  fl] 

Here  r  is  a  measure  of  instantaneous  risk  that  is  integrated  over  the  course  of  a  dive  and 
post-dive  period.  When  more  than  one  hypothetical  "tissue"  is  assumed  to  exist,  ws 
obtain  the  total  instantaneous  r  by  summing  the  col  hibutions  from  the  individual  tissues: 

r  ^  [2] 

i 

r  here >1  «  instantancaus  risk  due  to  tissue  i.  The  form  we  give  to  t,  !or  tissue  i  is  as 
follows: 


h  -  1  [  -  P«.b  -  Pthr  +  lco2(Po2)  ]  /  Pamb  .  fj  >  0  [3] 


where  Pjjjj  »  partial  pressure  of  inert  gases  in  tissue  i  (fsw); 

Pwnb  »  ambient  hydrostatic  pressure  (fsw); 

Po2  =  partial  pressure  of  oxygen  in  breathing  gas  (fsw); 

A,  =  gain  factor  for  tissue  i  (min’); 

Pjjj,  ■  threshold  pressure  difference  (fsw); 

ko2  “  a  risk  coefficient  for  Po2  (dimensionless). 

Thus,  when  two  or  more  "tissues"  are  postulated,  the  models  predict  that  the  probability 
that  DCS  will  occur  is  the  joint  probability  of  DCS  in  all  of  ±e  tissues.  The  metabolic 
gases  CO2  and  HjO  are  ignored  in  this  calculation.  Whenever  the  numerator  on  the 
right  side  of  [3]  is  less  than  zero,  r|  is  set  equal  to  zero,  so  that  the  integrated  risk  caimot 
diminish  with  time.  In  other  words,  risk  can  accumulate  but  it  carmot  be  depleted.  P^^, 
is  an  al  'olutely  safe  excess  partial  pressure  of  inert  gas  that  can  be  sustained  indefinitely 
with  no  risk  .^f  DCS.  A  non<zero  value  of  koj  J)dicates  that  the  risk  of  DCS  depends  on 
the  partial  pressure  of  0]  as  well  as  the  partial  pressure  of  the  inert  gas: 
kQ2  >  0  suggests  that  breathing  a  high  P02  increases  the  risk  of  DCS  independently  of 
the  effect  of  the  inert  gas,  and  ko2  <  0  suggests  that  a  high  P02  is  beneficial  from  the 
standpoint  of  DCS  prevention,  aside  firom  the  decreased  partial  pressure  of  inert  gas  that 
it  implies.  The  risk  model  contains  three  adjustable  parameters  besides  those  used  in 
computing  Ptis:  A,  koj ,  and  Pthr-  These  last  two  parameters  can  easily  have  values  of 
zero  and  can  be  fixed  at  zero  m  order  to  simplify  the  fitting  routine.  Tlie  parameter  A 
would  be  zero  only  if  there  was  no  risk  of  DCS,  regardless  of  the  dive  profile.  P,i^  in 
equation  [3]  is  calculated  by  assuming  that  gas  exchange  kinetics  in  the  hypothetical 
tissue  are  either  mono-expcnentid  or  bi-exponential.  The  various  models  consist  of 


combinations  ot  various  numbers  of  "tissues”  governed  by  one  or  the  other  of  these 
kinetic  models.  The  following  is  a  summary  of  the  models;  they  are  presented  in 
increasing  order  of  complexity. 

The  first  and  simplest  model  incorporating  gas  exchange  is  summarized  as  follows: 
Gas  exchange  model  1: 

single  mono-exponential  tissue,  having  time  constant  t ; 

1  parameter  in  gas  exchange  model:  t; 

1- 3  other  parameters:  A,  Pu.r  (?).  ko2  (?)• 

Thus,  the  risk  model  contains  a  total  of  2  to  4  parameters,  depending  upon  whether  Pth, 
and/or  Icqj  are  rsed.  A  logical  extension  of  the  single  mono-exponential  tissue  model  is 
a  model  including  two  such  tissues  in  parallel: 

Gas  exchange  model  2: 

two  mono-exponential  tissues,  having  tme  constants  and  Xj; 

2  parameters  in  gas  exchange  model:  x^  and  Xj; 

2- 4  other  parameters:  A  ^ ,  A  j ,  P®,  P)  koj  P). 

In  the  double-exponential  description  of  gas  exchange  in  a  single  tissue,  the  single 
exponential  is  replaced  by  the  sum  of  two  exponentials.  The  first  exponential  is 
multiplied  by  the  dimensionless  normalized  weighrmg  factor  W| ,  and  the  second 
exponential  is  multiplied  by  (l*Wj).  This  kinetic  model  has  3  kinetic  parameters  rather 
than  the  1  of  a  single  exponential.  If  only  one  such  tissue  is  postulated,  then  the 
following  model  results: 
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one  do”ble-exponential  tissue,  having  time  constants  Tj  and 
Wi  is  the  weighting  given  to  the  first  exponential  and  (l-Wj) 
is  the  weighting  on  the  second  exponential; 

3  parameters  in  gas  exchange  model:  x  j ,  ,  w^; 

1-3  other  parameters:  A  ,  (?).  ko2  (?)• 

If  we  postulate  three  single-exponential  tissues,  then  this  model  results: 

Gas  exchange  model  4: 

three  mono-exponential  tissues,  having  time  constants  xi ,  xj ,  and  Xj; 

3  parameters  in  gas  exchange  model:  x ^ ,  and 
3-5  other  parameters:  A  j ,  A  2 ,  A  3 ,  Pu.  (?).  ko2  (?). 

Two  double-exponmtial  tissues  result  in  the  following  model: 

Gas  exchange  model  5: 

two  double-exponential  tissues  A  and  B,  having  time  constants  x^i  >  x^  , 

and  X  02; 

w^i  and  Wgj  are  the  weighings  given  to  the  first  exponentials  in  tissues  A 
and  B;  (l-w^i)  and  (l-Wgj)  are  the  weighings  on  the  sectmd  exponentials  in 
tissues  A  and  B. 


6  parameters  in  gas  exchange  model:  ,  x^ ,  w^j ,  x^  ,  x^ ,  Wgi; 

2-4  other  parameters:  A  ,  A  g ,  (?),  ko2  (?)• 

The  means  of  computing  Ptis  are  detailed  in  Report  I  of  this  series  (1).  Note  that 
although  reference  has  been  made  to  recognizable  pl^ical  processes  in  formulating 


these  models,  the  models  are  not  to  be  taken  as  literal  representations  of  the  truth  and 
the  parameter  values  obtained  by  fitting  them  to  data  are  not  to  be  regarded  as  having 
significance  outside  of  their  use  with  these  models. 
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III.  DATA  SOURCES  AND  HANDUl^G 


Sumicaries  of  the  data  used  in  this  stud>  are  offered  in  Table  1.  Somewhat  more 
detailed  summaries  will  be  offered  in  a  future  report  The  dives  in  data  sets  whose 
names  begin  with  "DC  or  DD"  were  performed  at  the  Defense  and  Civil  Institute  of 
Experimental  Medicine  under  the  direction  of  R.Y.  Nishi  (6-11).  All  data  sets  whose 
names  begin  with  "EDIT  are  based  on  dives  at  the  Experimental  Diving  Unit  in  Tanama 
City,  Florida,  and  were  collected  imder  the  direction  of  EJD.  Thalmann  (12-14).  Data 
set  NMR8697  was  collected  at  the  Naval  Medical  Reseturh  Institute  under  the  direction 
of  PJC  Weathersby  (15).  All  of  the  above  data  are  from  wet  chamber,  working  dives. 
Diy  dives  were  excluded  because  of  the  possibility  that  immersion  is  one  of  the  factors 
that  controls  the  risk  of  DCS  (16).  In  almost  all  cases,  moderate  physical  work  done 
at  depth,  not  during  decompression. 

Automated  recordings  of  the  depth  'r.  ^e  profiles  were  available  for  all  of  the 
dives.  These  were  converted  to  the  format  suitable  for  our  analysis  by  use  of  a  computer 
algorithm  that  simplifies  each  depth/time  plot  to  a  sequence  of  up  to  76  connected  line 
segments.  This  simplified  profile  was  required  to  agree  with  the  original  depth/time 
recording  to  within  1  ft  in  depth  and  0.1  min  in  tirae;  also,  a  portion  of  the  original 
recording  was  considered  eligible  fo;  representatic.^.  as  a  single  line  segment  only  if  it 
was  linear  to  within  .0-15%.  Depth  of  water  in  a  suspended  wet  pot  was  taken  to  be  the 
height  of  water  above  mid-chest  level  of  a  diver  of  average  hei^t.  Divers  used  the 
Mark  15  or  Mark  16  breathing  apparatus  during  the  0.7  dives,  so  that  the  co^rjposition  of 


the  diver's  breathing  gas  was  regulated  via  the  continuous-feedback  control  afforded  by 
this  system,  which  typically  is  assumed  to  maintain  the  partial  pressure  of  Oj  within 
0.1  atm  of  the  set  point  The  divers  breathed  O2  by  mask  (9-11)  during  decompression 
in  the  DC8AOW  data  set  (the  gas  composition  was  assumed  to  be  99.5%  O2 , 0.5%  N2). 
Gas  changeovers  were  encoded  by  representing  the  gas  composition  as  a  linear  function 
of  time.  During  the  changeovers  between  O2  and  air  in  the  DC8AOW  data,  it  was 
assumed  that  the  simultaneous  washout  of  the  breathing  apparatus  and  the  lungs 
required  a  total  of  13  min;  all  other  gas  changeovers  in  all  other  data  sets  were  assumed 
to  require  1  min. 

A  post-dive  siuface  time  of  either  12  or  24  h  was  included  with  each  dive  profile. 
This  number  defines  the  upper  integration  limit  for  computing  the  risk  integrals  defined 
in  Part  IL  It  represents  the  minimum  length  of  time  spent  on  the  surface  between  dives, 
based  on  the  best  information  provided  by  the  people  who  directed  the  dives.  It  is  highly 
probable  that  the  tissue  supersaturadon  (and  hence  the  risk  function)  decreases  to  zero 
for  all  of  these  dives  before  12  h  has  elapsed  post-dive,  even  when  much  longer-than- 
risual  time  constants  are  used  in  the  models.  Questions  about  the  data  w  '.re  resolved 
through  consultation  with  the  men  in  charge  of  the  dives.  Each  dive  outcome  was  set  to 
zero  for  a  definitely  safe  dive,  03  for  marginal  symptoms  (skin  itch  or  discoloration,  mild 
pain  of  brief  duration,  and  moderate  fatigue),  and  1.0  for  any  more  severe  symptoms 
ascribed  to  DCS.  When  a  dive  outcome  was  seriously  in  doubt,  the  dive  was  excluded 
altogether  from  the  analysis.  Cases  of  recompression  without  symptoms  were  deemed 
cases  of  unknown  outcome  and  were  excluded,  becaus*.'.  it  was  considered  that  symptoms 
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could  have  appeared  if  the  dives  had  gone  to  completion.  For  those  cases  in  which  a 
stricken  diver  was  recompressed,  his  dive  was  encoded  as  though  the  intended  schedule 
had  been  completed.  This  convention  was  followed  out  of  necessity:  as  a  consequence  of 
the  use  of  cumulative  risk  mtegrals,  a  dive  that  had  been  truncated  for  recompression 
therapy  would  always  be  computed  by  any  of  our  models  as  being  safer  than  a  completed 
dive  on  the  same  schedule.  Thus,  since  the  truncated  dive  is  associated  with  a  case  of 
DCS  and  the  completed  dive  generally  is  not;  the  model  would  be  forced  to  fit 
dose-response  data  in  which  decreased  dose  apparently  results  in  high^ '  response.  We 
have  observed  that  the  fit  of  our  models  to  data  is  markedly  degraded  by  the  inclusion  of 
dive^  which  have  been  encoded  as  truncated  and  have  outcome  equal  to  1.  Although  the 
procedure  we  have  followed  seems  contrived,  it  is  the  most  rational  we  have  yet  devised. 

There  are  some  instances  of  repeat  dive  trials  in  whidh  a  diver  suffered  DCS  during 
the  surface  interval  or  before  completing  the  first  dive,  and  was  therefore  recompressed 
for  treatment  without  descending  for  the  second  dive  in  the  series.  These  were  encoded 
as  single  dives  and  were  put  into  data  sets  with  the  other  single  dives.  This  seems  to  us 
the  most  logical  way  to  treat  these  cases,  but  note  that  it  biases  the  data  by  making  the 
single  dive  sets  appear  more  hazardous  while  simultaneously  making  the  repeat  dive  sets 
tq>pear  safer.  Our  decision  about  how  to  organize  the  data  into  sets  changes  nothing 
when  we  fit  models  to  combined  single  and  repeat  data;  the  data  were  not  changed,  only 
their  arrangement  into  subsets  was. 
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IV.  DATA  ANALYSIS 


The  evaluation  of  the  data  in  Table  1  begins  with  their  being  fitted  individually  by 
the  models  described  in  Part  n  using  likelihood  maximization  (17).  Simply  put,  an 
iterative  nmnerical  algorithm  [more  specifically,  a  modified  Marquardt  nonlinear 
least-squares  minimization  algorithm  (18)]  is  employed  to  determine  the  set  of  parameter 
values  that  maximizes  the  probability  of  the  observed  family  of  outcomes.  Since,  in 
general,  the  global  maximum  will  be  surrounded  a  host  of  lesser  maxima  on  the 
likelihood  surface,  it  is  always  necessary  to  make  multiple  attempts  at  fitting,  using  a 
variety  of  initial  guesses  of  the  parameter  values,  before  one  can  be  reasonably  assured 
of  having  achieved  convergence  at  the  global  maximum.  Following  the  fitting  of 
individual  data  sets,  various  relevant  combinations  of  those  sets  were  also  fitted  by 
likelihood  maximization.  The  fits  to  these  combinations  were  compared  with  the  fits  of 
the  original,  smaller  data  sets,  and  this  allowed  us  to  apply  one  of  the  statistical  tests  for 
comparing  single  dives  with  repeat  dives  or  air  dives  with  0.7  dives,  lire  test  is  as 
follows  (17). 

Let  a  data  set  be  fitted  using  likelihood  maximization  by  two  models  such  that  the 
first  model  is  a  subset  of  the  second,  that  is,  the  second  model  consists  of  the  first  model 
plus  one  or  more  additional  adjustable  parameters.  Then,  if  the  second,  more  elaborate, 
model  is  not  intrinsically  superior  to  the  first  model  for  this  data  set,  then  the  test 
statistic  2(LL  2  •  UL  j)is  distributed  approximately  as 
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X^  V  =  (nj  -  ni)  [4] 

where  LL  ^  =  ln(ma;dmum  likelihood)  by  the  first  model; 

LL  2  =  ln(maximum  likelihood)  by  the  second  model; 

=  the  chi'Square  function  having  v  degrees  of  freedom; 

%  =  number  of  parameters  in  the  first  model; 
nj  »  number  of  parameters  in  the  second  model 

In  words,  if  the  two  models  are  equivalent  for  this  data  set,  then  twice  the  difference. 

in  tne  log  likelihoods  is  a  random  variable  approximately  distributed  as  the  chi  square 

probability  density  function  having  a  number  of  degrees  of  freedom  equal  to  the 

difference  in  the  munber  of  parameters.  Thus,  to  compare  the  models  we  reject  the  null 

hypothesis  that  they  are  equivalent  if  2(LL  j  *  LL  j)  U  an  improbably  large  variable 

with  V  a  (n2  -  nj)  degrees  of  freedom.  For  example,  if  (nj  •  nj)  *  3,  then  even  if  the 

models  are  equivalent  there  is  still  a  0.0S  probability  that  2(LL  2  >  LL  j)  will  be  greater 

than  5.99,  because  the  area  tmder  the  x^  curve  to  the  right  of  5.99  is  equal  to  0.05.  The 

above  procedure  is  the  likelihood  ratio"  test 

Now,  if  the  two  data  sets  are  each  fitted  by  a  3-parameter  model  and  then  their 

combination  is  also  fitted  by  the  same  model  we  can  consider  the  results  as  basis  for 

comparing  a  3-parameter  model  with  a  6-parameter  model  for  the  superset.  If  the  6- 

parameter  model  is  not  demonstrated  to  be  significantly  superior  to  the  3-parameter 

model  then  we  accept  that  the  two  subsets  are  from  the  same  population  and  call  them 

"combinable”.  The  sense  to  this  is  that  if  two  data  sets  can  be  "combined",  then  one 

model  can  describe  all  of  the  data  about  as  well  as  it  can  describe  the  individual  data 

sets;  This  is  a  measure  of  the  similarity  of  the  two  data  sets,  or  at  least  of  how  similar 

they  seem  to  the  model 
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A  second  statistical  test  for  similarity  between  two  data  sets  is  accomplished  by  first 
fitting  a  model  to  a  data  set  to  determine  the  optimal  parameter  values,  and  then  using 
the  same  model  with  these  parameter  values  to  predict  the  P(DCS)  of  the  dives  in  the 
other  data  set  The  predicted  average  P(IX^)  is  then  compared  with  die  observed 
average  P(DCS),  which  is  simply  equal  to  the  raw  incidence  of  DCS. 

The  uncertainty  on  the  predicted  P(DCS)  is  estimated  using  a  propagation  of  error 
that  follows  reference  #19.  Let  the  predicted  P(DCS)  be  represented  by  P\  then  P  can 
be  written  as  a  function  of  the  adjustable  parameters  in  the  risk  model  used  to  calculate 
P(DCS): 

P  =  P  (B, ,  i  =  1, 2,  3  ...Nprm)  [5] 

where  B,  »  an  adjustable  parameter, 

NpKM  *  total  number  of  adjustable  parameters  in  the  risk  model. 

Then  the  variance  on  the  estimate  of  P(DCS)  for  dive  k  is  determined  from  the 
covariance  matrix  as  follows: 

NpRM  Npmyj 

_  _  aPfc  aPk 

var  (PJ  »  Zrf  Z/  —  —  cov  B|  cov  Bi  [6] 

i«l  j«l  aB,  aBj 

The  average  predicted  P(tXrS)  for  an  entire  data  set  simply  equals  the  sum  of  P  for  all 
dives  divided  the  number  of  dives: 
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where 


[7] 


NdiV  » 


Npiv 

k=l 


where  Nqiv  =  number  of  dives  in  the  data  set  To  compute  the  variance  on  we 
write  an  expression  analogous  to  equation  [6]: 


Npiu4  NpRM 


var 


-  r  E 


cov  cov 


i»l  j=l  3B|  3Bj 
into  which  we  make  the  following  substitutions: 


Njjiv 

T, 

aF, 

dB, 

k=l 

dB, 

Npiv 

aF, 

dBj 

k»l 

aBj 

Then  from  equation  [7]  we  have 

var(^w,)  •  var  (F^)  /  (Nnn^)* 


[8] 


[9a] 


[9b] 


[10] 


When  converting  var(F^  to  a  confidence  interval  on  ,  we  assume  to  be 
normally  distributed. 


V.  DISCUSSION  OF  RESULTS 

A.  Comparing  and  Combining  Data  Sets 

Likelihood  maximization  was  used  to  St  models  to  the  original  data  sets  summarized 
in  Table  1  and  to  various  supersets  formed  by  combining  some  of  the  original  data  sets. 
Tables  2*7  contain  selected  restilts;  these  tables  list  the  In(likelihoods)  and  optimal 
parameter  values  obtained  via  the  above  procedure  for  certain  key  supersets. 

Appendices  1*10  contain  the  same  data-Stting  informaSon  for  all  of  the  sets  and 
supenets  to  which  likelihood  maximization  was  implied. 

The  aim  of  this  study  is  to  determine  whether  a  single  model  can,  with  satisfactory 
accuracy,  predict  the  decompression  risks  of  single,  repeat,  and  multilevel  dives  on  any 
N2/O2  breathing  mixture.  An  examination  of  the  table  entries  allows  one  to  draw 
several  conclusions  central  to  this  aim.  However,  the  useful  information  can  be  extracted 
only  by  comparing  entries  between  tables  as  well  as  within  each  table,  tbe  authors  well 
recognize  how  tedious  this  is  and  will  try  to  select  the  most  significant  results  for 
mention  in  the  main  text 

One  note  on  the  tabulated  results:  in  Tables  2*10  and  Appendices  1*10,  the  symbol 
7///"  indicates  that  the  parameter  values  converge  such  that  the  model  becomes 
equivalent  to  model  1,  the  single-tissue  monoexponential  model,  which  is  our  simplest 
risk  accumulation  model.  AU  of  the  more  complicated  risk  accumulation  models  reduce 
to  model  1  when  certain  of  their  parameters  equal  zero. 
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Table  2  lists  the  highest  In(likelihood)  values  obtained  for  the  largest  data  sets  that 
can  be  formed  without  mixing  data  from  different  "categories".  In  other  words,  single 
dive  data  are  not  mixed  with  repeat  dive  data,  dives  on  different  breathing  mixes  are  not 
mbced,  etc.  Data  set  EDU1180R,  which  contains  our  only  multilevel  dive  data,  has  been 
emitted  from  Table  2  and  from  all  subsequent  tables  oecause  it  is  not  "combinable"  witli 
other  data  sets  according  to  the  likelihood  ratio  test  (see  "Data  Analysis"  se^nion).  We 
will  elaborate  on  this  point  in  Part  B  of  this  section.  Table  3  lists  the  LL  values 
obtained  for  the  largest  possible  groupings  of  single  dive  data,  repeat  dive  data,  air  dives, 
and  0.7  dives,  where  data  from  different  categories  have  been  mixed.  Table  3  also 
contains  LL  values  obtained  for  still  larger  groupings  of  dives  from  all  categories. 

y^pendix  1  contains  the  In(likelihood)  values  obtained  for  all  data  groupings  in 
which  data  from  different  categories  has  not  been  mixed;  ^pendix  2  lists  the  LL  values 
for  combinations  of  data  from  different  Categories  and  Appendix  3  is  for  groupings  of 
data  from  all  categories.  The  results  in  Aprendix  1-3  were  obtained  speciffcally  to  assess 
the  combinability  of  various  groupings  of  aata  using  the  likelihood  ratio  test.  Only  the 
relatively  simple  models,  having  ho  more  than  4  adjustable  parameters,  were  applied 
becaiise  many  of  the  data  sets  are  so  small  that  there  is  no  advantage  to  using  more 
complex  models.  Note  (in  Appendix  1)  that  the  original  data  s.ts  EDUI180R, 
EDU88SM,  DDREPWET,  and  DC8AOW  were  not  dealt  with  individually  because  they 
had  too  few  cases  of  DCS  (no  more  than  4)  to  be  fitted  meaningfully  by  any  of  the 
models.  Strictly  speaking,  this  makes  it  impossible  to  assess  their  combinability  with 
other  data  sets  using  the  likelihood  ratio  test 


One  conclusion  is  possible  inunediatel/:  given  a  data  set  large  enough  to  support  the 
use  of  models  having  multiple  parameters,  it  is  always  true  that  at  least  one  of  the 
mo^jls  based  on  integrated  risk  fits  the  data  substantially  better  than  the  null  model.  As 
an  illustration,  note  the  entries  in  Table  2  for  single  air  dives.  For  the  null  model,  which 
has  1  adjustable  parameter,  the  ln(likeliho'  i)  is  -154.9.  Model  1  (the  1 
mono-exponential  model)  without  a  threshold  overpressure  is  a  2-parameter  model;  the 
maximum  LL  obtained  with  this  model  is  -153.9.  To  evaluate  tne  improvement  in  fit 
afforded  by  model  1,  we  make  use  of  the  likelihood  ratio  test  (equation  [4]),  which  uses 
the  chi-square  function  with  1  degree  of  freedom.  We  find  that  we  can  be  only  80% 
confident  in  considering  model  1  to  be  intrinsically  superior  to  the  null  model  for  these 
data.  Adding  the  threshold  overpresstire  to  mo'.el  1,  for  a  total  of  3  adjustable 
parameters,  fails  to  improve  the  fit  at  all.  However,  the  4-parameter  models  are 
considerably  more  effective:  model  2  (the  2  mono-exponential  model)  improves  on  the 
null  model  by  9.6  LL  units,  add  model  3  (with  1  bi-exponential)  is  7*2  LL  units  "better” 
than  the  null  model.  Both  results  are  significant  at  the  p<  0.005  level. 

Note  that  a  difference  in  LL's  means  that  two  models  differ  in  their  predictions  of 
the  likelihood  of  the  observed  family  of  dive  outcomes.  For  a  difference  of  9.6  LL  units 
this  difference  in  predicted  likelihoods  is  a  factor  of  exp(9.6)  >  15,300.  Since  the  hull 
model  assumes  a  uniform  risk  of  DCS  regardless  of  dive  profile,  the  superior  fit  of  our 
more  sophisticated  models  reflects  their  ability  to  discriminate  among  dives  with  regard 
to  their  riskiness. 
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For  the  data  listed  in  Table  2,  with  the  exception  of  the  single  air  dives,  noted 
above,  neither  of  the  4-parameter  models  provides  an  improvement  of  more  than  1.6  LL 
units  (p<025)  over  the  fit  cf  model  1.  For  the  generally  larger  and  more  diverse  data 
sets  in  Table  3,  we  see  that  model  2  usually  provides  the  best  fit.  For  example,  for  the 
most  complex  duta  s'.t,  containing  1878  single  and  repeat  dives  on  various  gas  mixtures, 
the  fit  of  model  2  is  11.7  LL  units  better  than  any  other  model's.  In  fact  for  the  data 
sets  in  Tables  2  and  3  we  find  that  model  2  always  provides  the  best  fit  for  data  sets 
having  more  than  500  dives,  but  never  improves  significantly  on  model  1  when  n<S00. 
This  illustrates  that  the  amount  of  information  one  gets  from  a  model  depends  on  how 
much  information  is  contained  in  the  data  set  to  which  it  has  been  fit;  a  simple  data  set 
cannot  support  an  elaborate  model. 

Tables  4<7  list  the  optimal  parameter  values,  with  their  standard  errors,  obtained  for 
the  selected  data  groupings  that  appear  in  Tables  2  and  3.  Tables  4  and  S  list  the 
optimal  parameter  values  obtained  by  fitting  the  1  mono-exponential  tissue  models, 
either  with  the  threshold  pressure  difference  treated  as  an  adjustable  parameter  (3 
adjustable  parameters)  or  with  P,hr  fixed  at  7.ero  (2  adjustable  parametesrs).  Tables  6  and 
7  contain  the  optimal  parameter  values  obtiuned  by  fitting  the  2  mono-exponential 
tissues  model  and  the  1  bi-exponential  tissue  model,  both  with  fixed  equal  to  '  ro  (4 
adjustable  parameters). 

Appendices  4-6  give  the  complete  listing  of  the  optimal  parameter  values  obtai'  ed 
for  the  2-  and  3-parameter  models.  Appendices  7-9  have  the  parameter  values  for  *he 
4-parameter  models. 
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An  important  observation  is  that  the  optimal  time  constant(s)  for  a  given  data  set 
depend  on  the  time  scale  of  the  dives  in  that  set  The  air  dives  (EDU885A, 
EDU885AR,  DC4W,  and  DRREPWET)  are  good  examples.  In  DC4W  the  dives  are  as 
short  as  0.5  h  and  never  run  as  long  as  2  h,  from  initial  descent  to  siufacing.  In 
EDU885AR  the  shortest  dives  last  almost  2  h,  most  run  over  3  h  and  more  than 
one-sixth  of  them  last  about  8  h.  EDU885A  contains  a  veiy  diverse  collection  of  dives 
lasting  anywhere  from  03  to  almost  7  h,  with  an  average  duration  intermediate  between 
those  of  the  other  two  data  sets.  The  significance  is  that  the  fitting  algorithm  cannot  use 
large  time  constants  to  much  advantage  when  describing  short  dives,  and  likewise,  small 
time  constants  will  not  be  used  in  fitting  long  dives.  (This  is  eqiiivalent  to  saying  that 
very  "slow"  tissues  are  unimportant  in  short  dives  and  that  "fast"  tissues  tend  not  to  affect 
the  outcome  of  long  dives,  if  one  assumes  each  time  constant  to  be  truly  associated  with 
a  distinct  tissue).  Again,  the  parameter  values  must  reflect  the  sort  of  information 
contained  in  the  data  for  which  they  were  optimized.  Not  surprisingly,  then.  Appendices 
4  and  7  show  that  DC4W  wants  the  shortest  time  cor,stants  and  EOU885AR  wants  the 
longest  ones.  Fcr  example,  for  the  single  tissue  mono-exponential  model  (model  1)  x  s 
46  min  for  DC4W,  t  »  281  min  for  EDU88SAR,  and  t  ■  173  min  for  EDU885A.  In 
view  of  the  disparity  in  fitted  parameter  values  among  the  air  data  sets,  it  is  not 
surprising  also  that  they  generally  cannot  be  combined  with  one  another  under  models 
having  just  one  or  two  adjustable  time  constants,  regardless  of  whether  we  are  comparing 
single  dive  data  sets  with  other  single  dive  sets  or  singles  with  repeats.  For  example,  by 
applying  the  likelihood  ratio  test  to  data  in  Appendices  1  and  2,  we  find  that  DC4W  and 
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EDU885A  (both  single  air  sets)  are  not  combinable  under  any  of  the  models,  that  data 
sets  EDU885A  (single  air)  and  EDU885AR  (repeat  air)  are  not  combinable  under  either 
model  2  or  model  3,  and  that  DC4W  does  not  combine  with  E33U88jA+EDU885AR 
according  to  any  of  the  models.  More  elaborate  models,  having  3  or  more  adjustable 
time  constants,  probably  could  accommodate  all  of  the  air  data  simultaneously,  but  no 
such  models  were  applied  to  the  air  data  exclusively. 

There  is  no  reason  to  conclude  that  there  is  any  fundamental  difference  in  the  dives 
from  one  air  dive  data  set  to  the  next,  other  than  their  duration.  The  optimal  parameter 
values  are  mathematical  devices,  not  precise  descriptors  of  definite  underlying  physical  or 
biological  processes. 

The  most  important  observation,  found  by  comparing  th6  entries  in  Tables  2  and  3, 
is  that  the  models  do  not  discriminate  between  single  and  repeat  dives  or  between  dives 
on  different  N2/O2  breathing  mixtures.  According  to  the  likelihood  ratio  test,  data  sets 
from  different  categories  generally  can  be  easOy  combined.  For  example,  using  model  2 
(the  2  mono-exponential  tissue  model,  4  adjustable  parameters),  which  appears  to  have 
the  widest  applicability  for  the  data  in  tips  smdy,  we  see  ttot  combining  data  sets  from 
two  categories  never  results  in  a  decrease  of  more  than  4.0  log  likelihood  (LL)  units  over 
the  sum  of  the  LL's  for  the  subsets.  Mathematically,  this  statement  can  be  written  (LL  2 
•  LL  1)  <  4.0,  tising  the  same  notation  as  in  equation  [4].  When  combining  single  air 
dives  with  repeat  air  dives  we  find  that  (LL  2  -  LL ,)  »  4.0 11.  units,  but  the  air  dive 
data  do  not  appear  to  combine  as  well  in  general  as  do  the  other  data,  as  was  discussed 
previously.  For  combining  single  0.7  dives  with  repeat  0.7  dives,  (LL  j  *  LL  ,)  =  0.5 


(keep  in  mind  that  the  data  set  EDU1180R  has  been  excluded  from  this  combination).  , 
For  the  combination  of  single  dives  on  both  air  and  0.7  (LL  2  -  LL  i)  =  1.1,  and  for  the 
combination  of  all  repeat  dives  (again  excluding  EDU1180R),  (LL  2  -  LL  j)  =  1.8.  By 
contrast,  the  likelihood  ratio  test  (equation  [4])  tells  us  that  if  two  data  sets  from  the 
same  population  are  described  u.sing  a  4-parameter  model,  there  is  still  a  0.05  probability 
that  the  highest  LL  that  will  be  obtained  for  their  combination  will  be  less  than  the  sum 
of  the  LL's  for  the  individual  sets  an  amount  of  at  least  4.74,  Le.,  (LL  2  -  LL  j)  > 

4.74.  Therefore,  the  LL  difference  must  be  at  least  4.74  before  we  would  reject  at  the 
0.05  level  the  null  hypothesis  that  single  and  repeat  dives,  or  air  and  0.7  dives,  are  from 
the  same  population.  Accordingly,  we  accept  that  the  likelihcod  ratio  test  does  nothing 
to  dissuade  us  from  the  nuU  hypothesis  that  single  and  repeat  dives  on  either  air  or  0.7 
atm  N2  are  described  equally  well  by  the  same  models. 

Table  9  summarizes  our  comparison  of  dives  from  different  categories  by  using  the 
models  as  predictors  'ather  than  descriptors,  as  outlined  in  the  Data  Analysis  section.  In 
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principle,  a  model  that  correctly  relates  P(DCS)  to  dive  proGles,  having  been  fitted  to 
some  data  set  ”A”,  should  be  able  to  accurat  ely  predict  the  outcomes  of  the  dives  in 
some  other  data  set  *B”.  The  realization  of  his  is  partly  confouixied  by  the  apparent 
randomness  of  dive  outcomes  (in  both  data  sets),  by  the  dependence  of  optimal  model 
parameter  values  on  the  specific  dive  profiles  contained  in  data  set  ”A^  and  by  the 
imperfections  of  the  models  themselves.  Nevertheless,  we  might  often  gain  insight  into 
how  similar  are  the  two  data  sets  by  comparing  the  predicted  average  P(DCS)  for  data 
set  ”B"  with  the  observed  average  P(DCS)  (vrhich  simply  equals  the  number  of  recorded 


DCS  cases  di\'ided  by  the  number  of  dives).  We  would  hope  to  find  substantial  overlap 
in  the  confidence  intervals.  The  confidence  intervals  on  the  PREDICFED  values  of 
average  P(DCS)  are  given  in  Table  9;  these  are  calculated  by  the  error  propagation 
method  shown  in  the  Data  Analy^s"  section. 

Unfortunately,  there  seems  to  be  no  satisfactory  way  to  calculate  a  confidence 
interval  on  the  estimate  of  average  P(DCS)  gr'en  by  the  OBSERVED  average  P(DCS). 
Many  of  the  dive  profiles  have  only  2  replicates  and  some  are  unreplicated,  and  binomial 
theory  canoOt  supply  useful  estimates  of  the  uncertainty  on  an  observed  proportion  when 
the  sample  size  is  only  1  or  2.  We  can,  however,  set  a  very  conservative  lower  bound  on 
this  uncertainty  by  assuming  that  all  of  the  dives  carry  the  same  underlying  P(DCS),  as 
though  the  data  set  is  composed  entirely  of  replicates  of  the  same  dive.  This  assumption 
leads  one  to  calculate  the  least  possible  uncertainly  in  the  estimate  of  P(DCS).  The 
confidence  interval  on  the  estimate  of  the  underlying  P(DCS)  by  the  observed  incidence 
is  easily  computed  usu^g  the  i^roximation  that  the  observed  P(DCJ)  is  a 
normally-Klistributed  random  variable. 

Note  that  in  Table  9a,  models  fitted  to  single  dive  data  are  used  to  predict  the 
outcomes  of  repeat  dives  and  models  fitted  to  air  dive  data  are  used  to  predict  the 
outcomes  of  0.7  dives,  whereas  in  Table  9b  the  situation  is  reversed.  Occasionally  the 
95%  confidence  interval  on  the  predicted  average  P(DCS)  does  miss  the  observed 
incidence,  as  for  example  when  models  fitted  to  the  repeat  0.7  dives  are  used  to  predict 
the  outcomes  of  repeat  air  dives  (see  Table  9b).  There  is  a  raw  incidence  of  72%  DCS 
in  the  repeat  air  data,  whereas  the  predicted  incidences  ranges  from  2.7  to  3.8% 


depending  on  the  mod^;l«  with  the  confidence  intervals  extending  no  higher  than  72%. 
However,  based  on  the  observed  incidence  o*  DCS  we  estimate  the  confidence  interval 
on  the  underlying  P(DCS)  to  span  £rom  3.6%  to  10.9%,  and  remember  xhat  this  is  a 
gross  underestimate  of  ±e  true  width  of  this  interval  In  summary,  we  recognize  no 
convincing  evidence  that  single  dives  are  incon^tibie  with  repeat  dives  or  that  air  dives 
and  0.7  dives  are  incompatible. 

B.  Multilevel  Dives 

Another  important  observation  is  that  the  addition  of  EDU1180R  to  a  data  set 
consistently  degrades  any  model's  ability  to  describe  the  data  in  that  set  EDU1180R 
consists  of  a  set  of  multilevel  profiles  in  which  the  diver  ascends  to  10-30  £sw  between 
descents.  Strictly  speaking,  one  caimot  assess  the  combinability  of  EDU1180R  with 
other  data  sets  without  fitting  the  models  to  EDU1180R  alone,  which  has  not  been  done 
because  of  the  sparseness  of  this  data  set  (2  DCS  cases  record^  out  of  128  man-dives). 
However,  if  one  applies  the  likelihood  ratio  test  on  the  assunqrtion  that  the 
risk-accumulation  models  would  fit  EDU1180R  jus*  as  well  as  the  null  model,  then 
EDU1180R  is  found  to  be  not  combinable  with  any  other  data.  Also,  the  addition  of 
EDU1180R  to  a  data  set  consistently  results  in  drastic  shortening  of  the  fitted  time 
constants  for  the  set 

Let  us  use  some  of  the  repeat  dive  data  as  an  illustration.  The  null  fit  to 
EDU1180R  yields  a  In(likelibood)  of  1030.  In  Appendix  2  we  see  that  the  combination 
of  data  sets  EDU184  and  EDU88SAR  can  be  fitted  with  model  2  so  that  UL  «  77.84 
(this  makes  them  very  easily  combinable  according  to  the  likelihood  ratio  test).  When 


EDU1180R  is  added,  the  best  fit  of  model  2  has  LL  »  9335.  This  much  degradation  of 
the  fit  is  significant  at  the  p<0.04  level.  In  i^pendix  8  we  see  that  the  optimized  time 
constants  for  [EDU184+EDU885AR]  are  t|  »  igg  niin  andTj  ^  min,  but  after 
EDU1180R  is  added  to  the  set  the  optimal  fit  is  obtained  with  «  0.43  min  and  = 
IDS  min.  We  get  similar  results  when  EDU1180R  is  mixed  with  other  data  sets. 

Table  8  summarizes  efforts  to  predict  P(DCS)  for  the  EDU1180R  dives  using 
parameter  values  optimized  for  other  data  sets.  Parameters  are  used  that  have  been 
optimized  for  nil  of  the  single  0.7  dives  (Le.,  EDU1180S  and  EDU88SM),  all  of  the 
repeat  0.7  dives  (mcludes  only  EDU184),  and  the  combination  of  all  dives  except 
EDU1180R.  In  no  instance  is  the  prediction  a  good  one,  with  between  11  and  32  DCS 
cases  predicted  as  compared  with  2  cases  observed.  The  poor  fit  is  also  indicated  by  the 
lew  LL  values  conqmed  with  the  null  model’s  LL,  indicating  that,  for  the  parameter 
values  obtained  by  fitting  other  data,  the  observed  dive  outcomes  are  far  fi-om  probable. 

The  explanation  is  that  EDU1180R  needs  considerably  shorter  time  constants  than 
do  the  other  data  sets  (with  the  exception  of  DC4W)  and  that  the  application  of 
longer-than-optimal  time  constants  is  unusually  deleterious  to  the  fit  Generally, 
lengthening  the  time  constants  in  one  of  the  models  has  the  partially  offsetting  effects  of 
reducing  the  gas  uptake  (thereby  reducing  the  total  risk  of  the  dive)  and  retarding  the 
offgassing  (thereby  inaeasing  the  total  risk).  This  dichotomy  acts  to  moderate  the  effect 
on  the  predicted  P(DCS)  of  changes  in  the  fitted  time  constants.  For  these  dives, 
however,  it  seems  that  the  former  effect  is  overwhelmed  by  the  latter,  so  that  there  is  a 
strong  intolerance  for  longer  time  constants.  Why  is  this?  The  dives  in  EDU1180R 
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contrast  with  our  other  data  in  that  the  divers  remain  underwater  for  4-S  h  before 
beginning  decompression  but  then  undergo  a  decompression  of  only  1-2  h  duration. 

Thus,  the  time  spent  imderwater  is  unusually  long  compared  with  the  decompression 
time,  so  that  even  "slow”  tissues  have  plen^  /  of  time  to  accumulate  inert  gas  but  only 
"fast”  tissues  can  offgas  effectively  before  surfacing.  In  our  models,  a  time  constant  much 
longer  than  about  60  min  results  in  the  prediction  of  a  substantial  tissue  supersaturation 
remaining  at  the  end  of  the  dive,  with  especially  severe  risk  accumulation  commencing  at 
the  20-  or  10-foot  stop  or  upon  surfacing  (depending  on  the  time  constant ).  Figure  1 
illustrates  this.  It  is  a  plot  of  one  of  the  multilevel  dives  in  which  ^**e  drers  descended 
to  150  fsw  and  stayed  there  for  29  min,  ascended  to  30  fsw  and  stayed  for  120  min,  and 
descended  again  to  150  fsw  for  30  min.  The  depth,  the  partial  pressure  of  N2  in  each  of 
the  "tissues”  postulated  by  model  2,  and  the  accumulated  P(DCS)  are  plotted  as  functions 
of  time.  The  partial  pressures  are  predicted  by  model  2  using  the  parameter  values 
optimized  for  aU  of  the  single-  and  repeat-  dive  data.  It  is  seen  that  the  tissue  having  the 
33-min  time  constant  (the  "fast”  tissue)  really  does  expel  most  of  its  dissolved  Nj  during 
the  stay  at  .30  ft,  and  consequently  this  tissue  makes  a  relatively  small  contribution  to  the 
P(DCS).  ^  contrast,  the  tissue  with  the  715-min  time  constant  actually  experiences  a 
small  net  gain  of  dissolved  during  the  30-£sw  stop,  and  the  final  decompression  is  far 
too  rtq>id  to  allow  a  thorough  washout  of  dissolved  gas  before  surfacing.  The  resulting 
overpressure  persists  for  6  h  versus  about  1  h  of  overpressure  in  the  fast  tissue.  In 
addition,  the  gain  coefficient  coupled  with  the  longer  time  constant  is  five  times  greater 
than  the  gain  coefficient  for  the  short  time  constant  (see  equation  [3]  for  the  definition 


of  the  gain  coefficient).  The  end  result  is  that  this  dive  profile  is  predicted  to  carry  a 
probability  of  DCS  with  the  715-min  tissue  accounting  for  the  bulk  of  the  risk. 

Similar  plots  could  be  shown  for  the  rest  of  the  multilevel  dives.  Therefore,  when 
EDU1180R  is  mixed  with  other  data  sets  and  the  combination  is  fitted  by  one  of  our 
models,  the  use  of  long  time  constants  results  in  a  gross  overestimation  of  P(DCS)  for 
the  EDUllSOR  dives,  and  the  use  of  short  time  constants  results  in  a  poor  fit  to  the 
remaining  data.  Of  course,  the  optimized  parameter  values  for  a  combination  of  data 
sets  never  quite  equal  the  optimal  values  for  any  one  of  the  subsets.  However,  the 
subsets  usually  appear  combinable  nonetheless,  using  "compromise”  parameter  values 
that  usually  are  within  uncertainty  bounds  of  the  parameter  estimates  for  the  individual 
subsets.  The  difference  with  EDUllSOR  is  the  unusual  sensitivity  of  the  fit  to  time 
constant  values,  and  this  sensitivity  is  a  consequence  of  the  special  characteristics  of 
these  dives. 

Whereas  the  fitted  parameter  values  are  not  taken  to  have  any  definite  physical 
significance,  ffie  only  justifiable  conclusion  is  that  a  fimdmental  difference  exists 
between  data  set  EDUllSOR  and  the  others  in  this  study  with  respect  to  their 
description  by  these  particular  models.  If  we  consider  the  models  to  be  suitable  for 
modelling  non-multilevel  dives,  then  it  is  necessary  to  doubt  whether  they  are  adequate 
for  modelling  multilevel  dives,  on  the  evidence  available  to  us  from  this  one  dive  series. 
C  Repetitive  Dive  Features 

An  evaluation  that  we  made  of  repeat  dives  may  be  of  interest  to  the  reader.  In  this 
exercise,  we  used  the  models  as  predictors  of  P(D^)  for  various  repeat  dive  profiles  and 
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varied  the  duration  of  the  surface  interval  as  an  independent  parameter,  without 
changing  any  of  the  other  aspects  of  the  dive  profiles.  We  found  that  P(DCS)  is  partly 
determined  by  two  effects  peculiar  to  repeat  dive  combinations.  First,  the  descent  for 
the  second  dive  eliminates  the  tissue  supersaturation  persisting  from  the  first  dive, 
thereby  halting  the  accumulation  of  risk  from  the  first  dive  and  reducing  P(DCS)  for  the 
dive  combination  as  a  whole.  This  is  effectively  a  truncation  of  the  first  dive.  Note  that 
this  reduction  in  P(DCS)  arises  from  our  assumption  of  an  accumulating  risk,  rather  than 
a  risk  resulting  from  some  single  instant  of  supersaturation  occurring  during  the  dive. 

The  second  effect  pecuh'ar  to  repeat  dives  is  the  addition  of  risk  to  the  second  dive 
because  of  the  residual  inert  gas  supersaturation  from  the  first  dive.  This  effect 
inaeases  the  risk  of  the  dive  combination. 

According  to  our  models,  for  the  repeat  dives  in  this  study,  the  lowest  P(DCS}  is 
achieved  with  a  surface  interval  of  infinite  duration,  i.e.  two  single  dives.  This  is 
another  way  of  saying  that  the  second  effect  described  above  dominates  the  first 
G>nsequently,  according  to  our  models  repeat  dives  are  inherently  more  dangerous  than 
single  dives  and  one  should  adjust  the  decompression  time  of  the  second  dive  to 
compensate.  As  an  example.  Figure  2  shows  the  predicted  P(DCS)  of  one  of  the 
double-dives  in  data  set  EDU8S5AR  as  a  function  of  surface  interval  time,  according  to 
model  2  with  its  parameter  values  optimized  for  all  of  the  repeat  dive  data.  The 
monotonic  decrease  with  siulace  interval  time  of  P(DC5)  is  typical,  although  maxima 
and  minima  are  possible.  Incidentally,  the  actual  dive  was  carried  out  with  an  86-minute 
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surface  interval;  there  was  a  98-minute  decompression  following  the  first  descent  and  a 
208-minute  decompression  after  the  second  descent 

It  certainly  is  possible  to  devise  dive  sequences  that  are  exceptions  to  the  rule,  that 
i&,  which  are  made  safer  by  a  reduction  in  the  times  allowed  between  descents.  One 
example  consists  of  a  hazardous  dive  followed  by  an  extremely  safe  dive  having  an 
extremely  long  decompression  time  -  this  is  equivalent  to  recompression  therapy. 
Another  example  occurs  when  the  second  dive  is  a  saturation  dive,  in  which  case  the 
residual  dissolved  inert  gas  from  the  first  dive  is  irrelevant. 


One  more  test  of  our  models  consists  of  attempting  to  improve  the  fit  by  allowing 
the  dimensionless  parameter  ko2  (see  equation  [3])  to  be  an  adjustable  variable  in  the 
curve-fitting  routine.  In  all  modelling  discussed  to  this  point,  Icqj  has  been  fixed  at  zero. 
If  the  fit  can  be  improved  significantly  by  allowing  to  be  non-zero,  it  would  suggest 
that  risk  accumulation  depends  on  the  partial  pressure  of  Oj  in  the  breathing  mix  as  well 
as  on  the  iartial  pressures  of  any  inert  gases.  Accordingly,  the  3-mono-exponential  tissue 


model  (model  4)  was  fitted  to  all  of  the  available  data  (except  EDU1180R)  with  k^j 
allowed  to  float  Table  10  and  Appendix  10  show  the  highest  Ih(likelihood)  and  optimal 
parameter  values  obtained,  and  these  can  be  compared  with  the  results  found  for  model 
4  with  kQ2  fixed  at  zero.  The  optimal  value  of  kos  is  scarcely  different  from  zero  (034 
with  a  standard  error  of  031)  but  the  improvement  in  fit  afforded  by  the  extra  parameter 
is  only  03  LL  units,  which  is  not  significant  We  conclude  that  we  have  no  strong 
evidence  to  support  the  inclusion  of  Fq]  in  our  risk  formulations.  Model  4  was  chosen 
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for  this  test  because  it  is  superior  to  our  other  models  at  fitting  all  of  the  data 
simultaneously. 
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VI.  MORE  COMPLEX  MODELS;  CONCLUSIONS 

The  siimmaiy  of  the  discussion  is  that  1)  our  models  have  some  substantial  face 
/alue  as  quantitative  descriptors  and  predictors  of  DCS  incidents,  2)  we  see  no 
compelling  reason  not  to  combine  data  regardless  of  whether  it  consists  of  single  dives  or 
repeat  dives  and  regardless  of  the  composition  of  the  Nj/Oj  breathing  mix,  and  3)  we 
see  no  compelling  evidence  that  any  of  the  data  sets  other  than  EDU1180R  comes  from 
a  different  population  than  the  bulk  of  the  data.  All  of  these  conclusions  are  based  on 
our  use  of  models  1,  2,  and  3.  It  remains  for  us  to  explore  whether  more  complicated 
functions  can  be  used  to  advantage.  Table  10  and  Appendix  10  each  show  the  highest 
In(likelihood),  and  the  optimized  parameter  values,  obtained  by  fitting  model  4  (the 
3*mono-exponential  tissue  model)  and  model  5  (the  2  bi-exponential  tissue  model),  each 
with  P«„=0  and  ko2^0,  to  all  of  the  available  data  (except  EDU1180R).  Model  4,  with 
6  adjustable  parameters  altogether,  actually  provides  a  somewhat  better  overall  fit  than 
does  the  8‘parameter  model  5.  The  excessively  high  standard  errors  on  the  model  5 
parameters  result  from  strong  cross-correlations  b^.tween  parameters,  suggesting  that  the 
data  set  may  not  be  sufficiently  large  and  diverse  to  support  a  model  as  elaborate  as  this. 
The  fit  of  model  4  is  seen  to  be  a  statistically  significant  improvement  over  that  of  model 
2  (the  2-mono-exponential  model,  having  4  parameters),  according  to  the  likelihood  ratio 
test. 

It  would  seem  that  the  3-mono-exponential  model  with  the  parameter  values  listed 
in  Table  10  represents  our  best  mathematical  description  of  the  available  data.  Its  fit 


improves  upon  the  null  model's  nt  by  20  LL  umts,  meamng  that  it  predicts  the  observed 
family  of  dive  outcomes  to  be  exp(20)  »  4.7*10^  times  more  probable  than  the  null 
model  does.  The  optimized  time  constants  of  0.^4  min,  129  min,  and  767  min  reflect  the 
diversity  of  the  available  data.  It  should  not  be  extrapolated  recklessly  to  dives  that  are 
appreciably  different  from  the  ones  in  our  database.  That  is,  it  should  not  be  expected 
to  be  a  reUabie  predictor  of  the  risks  of  dives  having  time  scales  different  from  those  in 
our  database;  as  we  have  seen,  the  optimal  time  constants  obtained  for  long  dives  do  a 
poor  job  of  describing  short  dives,  and  vice-versa.  Our  database  contains  dives  ranging 
from  03  to  8  h,  so  saturation  dives  are  not  encompassed.  Obviously,  it  should  not  be 
extrapolated  to  multi-level  dives,  which  were  excluded  from  the  data  set  for  which  the 
parameter  values  in  Table  10  were  optimized. 
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APPENDIX  1 


Log  Likelihoods  for  the  2-,  3-,  and  4-Paraaeter  Models: 

Date  Sets  Not  Combined  Across  Categories;  k(]2  ~  0 

In  the  following  tables,  "////"  indicates  that  the  model  converges  to  a 
solution  equivalent  to  a  two-parameter  DLE6  solution. 

-In  (likelihood) 

DCS  . - . 


cases/dives 

null 

model 

.  1 

model  2 

model  3 

single  air 

dives : 

no  Ptfc, 

W/Ptht 

EDU885A 

30/483 

112.41 

109.92 

//// 

101.88 

105.37 

DC4V 

18+4(0.5)1/244 

41.738 

40.595 

//// 

35.514 

37.054 

EDU88SA 
+  DC4U 

[38+4(0.5)1/727 

154.88 

153.86 

//// 

145.25 

147.70 

single  0.7 

dives : 

EDU1180S 

10/120 

34.420 

30.495 

30.474 

//// 

30.456 

EDU88SH 

4/81 

15.932 

.... 

.... 

.... 

.... 

EDU1180S  ■¥ 

EDU885M  14/201 

50.800 

44.851 

//// 

44.821 

//// 

single  dives,  other  gas  mixes: 

NMR8697 

[11+18(0.5)1/477 

83.010 

81.970 

81.939 

81.844 

81.931 

DC8A0W 

[2+1(0.5)1/45 

9.655 

.... 

.... 

.... 

repeat/multilevel  air  dives: 

EDU885AR 

11/182 

41.528 

38,027 

36.447 

//// 

36.390 

DRREPWET 

3/12 

6.748 

.... 

- - 

.... 

.... 

EDU88SAR  * 

DRREPWET  14/194 

50.285 

45.719 

44.240 

//// 

44.115 

repaat/mulcllevel  0.7  dives: 

EOU1180R 

2/li8 

10.302 

.... 

.... 

— 

.... 

EDU184 

11/234 

44.369 

38.455 

37.420 

37.930 

37.329 

EOU1180R  + 

EDU184  13/362 

56.011 

52.989 

52.914 

//// 

52.511 

APPEHDXZ  2 


Log  Likelihoods  for  the  2*.  3-,  end  4>Pereaeter  Models: 
&>abinetions  of  Data  Seta  Fron  Different  Categories 


-In  (likelihood) 


DCS 

caaes/dlves 

null 

model  1 

model  2 

model  ' 

no 

W/Ptht 

single,  air  -»■  0.7: 

EDU1180S  +  EDU885A 
*  EDU83SM  +  DC4V 

[52+4(0.5)1/928 

205.98 

204.17 

//// 

191.21 

195.46 

single,  any  gas  nix: 

EDU1180S  +  EDU885A 
+  EDU885M  +  DC4W 
+  DC8A0W  +  NMR8697 

(65+23(0.5)1/1450 

299.51 

297.23 

//// 

285.39 

//// 

repeat/eultilevel,  air  +0.7: 

EDU184  +  EDU885AR  22/416 

86.080 

78.298 

75.109 

77.837 

75.923 

EDU184  +  EDU885AR 
+  DRREPVET  25/428 

95.261 

86.869 

83.544 

86.442 

83.523 

EDU1180R  +  EDU184 

+  EDU885AR  24/544 

98.364 

93.459 

92.988 

93.346 

//// 

air,  single  +  repeat/nultllevel : 

E0U885A  +  EDU885AR  41/665 

153.94 

148.63 

//// 

144.19 

146.30 

EDU885A  +  EDU885AR 
+  DC4V  +  DRREPVET 

1 52+4(0. 5) 1/921 

205.55 

202.28 

//// 

194.95 

197.98 

EDU883A  +  EDU885AR 
+  DC4W  [49+4(0.5)1/909 

196.45 

194.23 

//// 

185.98 

189.43 

0.7,  single  +  repeet/aultilevel: 

EDU1180S  +  EDU1180R  12/248 

48.047 

47.500 

41.267 

46.656 

//// 

EDU1180S  +  EDU184  21/354 

EDU1180S 

+  EOU1180R  +  E0U184  23/482 

79.685 

92.419 

69.096 

87.992 

68.571 

86.962 

68.619 

//// 

68.287 

//// 

EDU1180S  + 

EDU184  +  E0U885M  25/435 

101.01 

83.583 

//// 

83.216 

83.395 

EDU1180S  +  EDU1180R 
+  EDU184  +  ZDU885M  27/563 

108.35 

104.89 

//// 

//// 

//// 
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APPENDIX  3 


Log  Likelihoods  for  the  2*,  3>.  and  4-Parameter  Models: 
Any  Category  of  Dive 


-In  (likelihood) 


DCS  . 

eases/dives  null  model  1  model  2 


any  type  of  dive: 


no  Ptfc,  W/Pthr 


EDU885A  -!• 

EDU885AR  +  EDU885M  45/746  169.98  165.00  ////  159.53 

EDU1180S 

+  ECU184  -I'  EDU885A 
E0U885AR  +  EDU885M 
*  0C4V  +  NMR8697 
+  DRREPtfET 

+  DC8A0W  (90+23(0.5) 1/1878  394.87  383.79  ////  377.09 


model  3 


162.04 


//// 


39 


APPENDIX  4 


Optialzed  Paraneter  Values  for  the  2>  and  3>Parameter  Models: 
Data  Sets  Not  Conbined  Across  Categorler 


Bodel  1 


(1  Tissue 
Mono - exponential ) 

Pthr  -  0 


model  1 
(1  Tissue 
Mono - exponential ) 
Pthr  0 


single  air  dives: 
EDU88SA 

EDU885A  -I-  DC4V 
DC4V 

EDU88SA  DC4W 

single  0.7  dives: 
EDU1180S 

EDU88SM 

EOU1180S  EDU88SM 


t  -  173  (38) 

A  -  0.00221  (0.00044) 

t  -  89.6  (16.8) 

A  -  0.00220  (0.000353) 

T  -  46.1  (28.2) 

A  -  0.00209  (0.00072) 

T  -  92.4  (16.8) 

A  -  0.00222  (0.00035) 


t  -  490  (140) 

A  -  0.00489  (0.00274) 


T  -  483  (88.6) 

A  -  0.00586  (0.00260) 


//// 


//// 


//// 


//// 


t  -  374  (516) 

A  -  0.00680  (0.0129) 

Pth,  -  2.08  (11.6) 


//// 


T  -  72.9  (28;7) 

A  -  0.00269  (0.00414) 
Ptt,  -  3.03  (11.2) 


t  -  134  (47.9) 

A  -  0.0816  (0.0978) 
Pth,  -  9.96  (3.98) 


t  -  177  (68.7) 

A  -  0.0491  (0.0703) 
Pth,  -  6.92  (4.35) 


single  dives,  other  gas  mixes: 

NMR8697  t  -  77.9  (34.2) 

A  -  0.00180  (0.000461) 

DC8A0W  . 

repeat/multllevel  air  dives: 

EDU885AR  T  -  281  (88.5) 

A  -  0.00406  (0.00249) 

DRREPWET  . 

EDU885AR  +  DRREPWET  t  -  297  (83.1) 

A  -  0.00510  (0.00280) 
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repeat/nultllevel  0.7  dives: 

EDU1180R  .  . 

EDU184  X  -  413  (90.0)  t  -  187  (66.5) 

A  -  0.00364  (0.00194)  A  -  0.0102  (0. 

Pthr  -  5.34  (3. 

EDU1180R  +  EDU184  t  -  101  (87.0)  t  -  104  (69.3) 

A  -  0.000875  (0.000256)  A  -  0.00118  (0 

-  1.92  (6. 


0118) 

37) 


.00163) 

96) 


APPENDIX  5 


Optimized  Parameter  Values  for  the  2-  and  3 -Parameter  Models: 
Combinations  of  Dives  From  Different  Categories 


model  1 
(1  Tl8st.e 

Mono-exponential) ; 
Pthr  -  0 


model  1 
(1  Tissue 
Mono - exponential ) 


'thr 


#  0 


single,  air  0.7: 

EOU1180S  +  EOU88SA 

+  EDU885M  +  DC4U  t  -  101  (15.1)  //// 

A  -  0.00223  (0.000307) 


single,  any  gas  mix: 

EOU1180S  +  EDU885A 

EDU885M  +  DC4W  x  -  99.6  (13.9)  //// 

•f  DC8A0U  -»■  NMR8697  A  -  0.00214  (0.000255) 


repeat/multilevel , 
E0U184  EOU885AR 


+  0,7: 

t  -  282  (67.0) 

A  -  0.00272  (0.00105) 


t  -  162  (29.0) 

A  -  0.0184  (0.0149) 
Pthr  -  6-6^  (2.03) 


EDU184  EDU885AR 
DRREPWET 


E0U1180R  +  EDU184  . 

*  EDU885AR 


T  -  288  (66.2) 

A  -  0.00307  (0.00115) 


T  -  93.9  (58.4) 

A  -  0.00105  (0.000216) 


t  -  160  (27.9) 

A  -  0.0218  (0.0165) 
Ptfc,  -  6.85  (1.90) 

T  -  92.6  (31.8) 

A  -  0.00279  (0.00321) 
Ptte  -  5.14  (4.G2) 


air,  single  4  repeat/multilevel: 
E0U885A 


4  FJ}U885AR 

t  -  186  (34,5) 

//// 

■ 

A  -  0.00228  (0.000436) 

EDU885A  4  DC4V 

t  -  103  (15.9) 

//// 

4  E0U885AR 

A  -  0.00202  (0.000293) 

EOU885A  4  EDU885AR 

t  -  103  (15.8) 

//// 

4  DC4V  DRREPVET 

A  -  0.00210  (0.000297) 
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0.7,  single  +  repeat/multilevel: 

EDU1180S  +  EDU1180R  t  -  126  (169) 

A  -  0.00112  (C. 000658) 


EDU1180S  -f  EOU184 


EDU1180S  + 
EDU1180R  -i-  E0U184 


EDU1180S  + 

EDU184  E0U88SM 


EDU1180S  +  EOU1180R 
■f  E0U184  +  EDU885M 


T  -  437  (65.2) 

A  -  0.00413  (0.00145) 

T  -  201  (100) 

A  -  0.00119  (0.000283) 

t  -  454  (52.9) 

A-  0.00486  (0.00153) 

T  -  231  (96.7) 

A  -  0.00135  (0.000431) 


T  -  73.7  (25.3) 

A  -  0.0575  (0.0799) 
Ptto  -  15.6  (1.60) 

t  -  232  (100) 

A  -  0.00700  (0.00488) 
Ptte  -  3.91  (3.22) 

T  -  172  (36,6) 

A  -  0.00740  '0.00497) 
Ptta  -  5.09  U.34) 


//// 


//// 


APPENDIX  6 


Optimized  PaTjameter  Values  for  the  2-  and 

Any  Category  of  Dive 

3. Parameter  Models; 

model  1 
(1  Tissue 

Mono -exponential) ; 

Pthr  -  0 

model  1 
(1  Tissue 

Mono - exponential ) 

Pthr  0 

EDU88SA  + 

EDU885AR  +  E0U88SM 

X  -  202  (31.7) 

A  -  0.00239  (0.000446) 

//// 

EDU1180S 

EDU184  +  EOU885A 
EDU88SAR  -I-  EDU883M 
+  DC4W  +  NMR8697 
■¥  DRREPVET  -t-  DC8A0V 

t  -  112  (12.7) 

A  -  0.00194  (0.000202) 

//// 

APPENDIX  7 


Optimized  Parameter  Values  for  the  A-Parameter  Models; 
Data  Sets  not  Combined  Aross  Categories 


model  2  model  3 

(2  Tissue  (1  Tissue 

Mono - exponential ) ;  B1 - exponential ) ; 

Pthr  -  0  Ptht  "  0 


single  air  dives: 


EDU88SA 


T,  -  22.5  (30) 
t,  -  730  (199) 

Ai  -  0.00173  (0.0016) 
A2  -  0.00940  (0.0075) 


12.6  (41.1) 

292  (44.6) 
0.822  (0.569) 
0.0296  (0.0237) 


DC4V 


t,  -  0.389  (0.149) 

Tj  -  358  (182) 

A,  -  0.553  (0.85) 

Aj  -  0.00970  (0.00014) 


T,  -  1.56  (4.2) 
t,  -  160  (158) 
w,  -  0.991  (0.031) 
A  -  0.0576  (0.023) 


EOU885A  +  DC4W 


-  27.3  (14.6) 

T,  -  749  (195) 

A,  -  0.00183  (0.00060) 
Ajt  -  0.00970  (0.00083) 


t,  -  14.5  (7.6) 

Tj  -  290  (58) 
w,  -  0.894  (  ) 

A  -  0.00921  (0.00026) 


single  0.7  dives: 
EOU1180S 

E0U885M 

EDU1180S  -»■  EDU885M 


//// 


•s,  -  2.51  (310) 

T,  -  493  (250) 

A,  -  0.00555  (1.07) 

Aj  -  0.00587  (0.00321) 


T,  -  12.7  (221) 
tj  -  375  (516) 

Wi  -  0.886  (1.66) 

A  -  0.00826  (0.0287) 


//// 


single  dives,  other  gas  mixes; 


NMR8697 


t,  -  0.263  (0.955) 

Tj  -  86.7  (54.1) 

A,  -  0.205  (2.57) 

A,  -  0.00166  (0.00120) 


T,  -  1.36  (31.2) 

T,  -  74.7  (60.7) 
vx  -  0.963  (0.718) 

A  -  0.00390  (0.0119) 


DC8A0V 


repeat/multllevel  air  dives: 

EDU885AR  ////  t,  -  120  (68.3) 

t,  -  2140  (3520) 
w,  -  0.892  (0.203) 
A  -  0.209  (0.339) 


DRREPVET 


EDU885AR  +  DRREPWET  ////  t,  -  141  (92.5) 

t,  -  2900  (7910) 
w,  -  0.925  (0.233) 
A  -  0.160  (0.264) 

repeat/ntultllevel  0.7  dives: 

EDU1180R  .  . 


EDU184 

t,  -  110  (634) 

T,  -  69.3  (54.3) 

T,  -  609  (282) 

t,  -  854  (334) 

Ai  -  0.000440  (0.00146) 

w,  -  0.733  (0.112) 

Aa  -  0.00642  (0.00515) 

A  -  0.0237  (0.0374) 

EDU1180R  +  EDU184 

//// 

T,  -  121  (117) 

T,  -  8820  (2.65E5) 

w,  -  0.991  (0.282) 

A  -  0.00429  (C. 00655) 


AFFENDZX  8 


Optimized  Parameter  Values  for  the  4>Parameter  Hc^els; 
Combinations  of  Data  from  Different  Categories 


model  2 
(2  Tissue 

Mono -exponential) ; 


‘thr 


"  0 


model  3 
(1  Tissue 
Bl- exponential) ; 


^thr 


a  0 


single,  air  4-0.7: 


EDUllSOS  4-  EDU885A 
4-  EDU885M  4-  DC4W 


T,  -  27.1  (12.6) 
t,  -  711  (147) 

A,  -  0.00183  (0.000583) 
A,  -  0.00872  (0.00539) 


t.  -  13.0  (5.35) 
tj  -  312  (39.5) 
w,  -  0.900  (0.0461) 

A  -  0.0102  (0.00249) 


single,  any  gas  mix: 

EDU1180S  4-  EDU88SA 

4-  EDU885M  4-  DC4W  t,  -  25.6  (10.3)  //// 

+  DC8AOW  +  NMR8697  t.  -  729  (115) 

A,  -  0.00204  (0.000534) 

Aa  -  0.00846  (0.00462) 


repeat/multilevel,  air  4-0.7: 


EDU184  4-  EDU885AR 


EDU184  4-  EDU885AR 
4-  DRREPWET 


T,  -  198  (270) 
tj  -  741  (782) 

A,  -  0.00143  (0.00248) 
A*  -  0.00525  (0.0107) 

T,  -  0.304  (0.109) 

T,  -  322  (90.4) 

A,  -  2.58  (10.8) 

Aj  -  0.00331  (0.00140) 


T,  -  148  (146) 
t,  -  2690  (11600) 
w,  -  0.932  (0.0224) 
A  -  0.0347  (0.0346) 

T,  -  0.755  (0.813) 

Tj  -  183  (40.1) 
w,  -0.995  (0.00570) 
A  -  0.0351  (0.0332) 


EDU1180R  4-  EDU184 
4-  E0U885AR 


T,  -  0.429  (1.01) 

Tj  -  105  (84.7) 

A,  -  0.108  (1.23) 

Aj  -  0.000992  (0.000274) 


//// 


air,  single  +  repeat/multllevel: 


EDU885A 
4-  EDU885AR 


t,  -  55.3  (51.4) 

T,  -  740  (217) 

a:  -  0.00105  (0.000271) 

Aj  -  0.00827  (0.00693) 


21.1  (49.2) 

351  (89.5) 

0.793  (0.410) 
0.00620  (0.00274) 


EDU88SA  4-  EDU885AR 
4-  DC4W  4-  DRREPWET 


t,  -  39.6  (20.2) 

Tj  -  732  (212) 

A,  -  0.00148  (0.00314) 
Aj  -  0.00807  (0.00677) 


t,  -  13.2  (7.28) 
t,  -  249  (63.6) 
w,  -  0.864  (0.0749) 

A  -  0.00624  (0.00191) 


EDU885A  4-  EDU885AR  r.  -  38.3  (19.8).  t,  -  12.5  (6.27) 

4-  DC4W  Tj  -  734  (200)  .  t,  -  258  (61.0) 

A,  -  0.00138  (0.000309)  w,  -  0.878  (0.0662) 
Aj-  0.00838  (0.00675)  A  -.0.00647  (0.00195) 


0.7,  single  +  repeat/miltllevel: 

,  . 

EDU1180S  +  EDU118bR 

//// 

//// 

EDU1180S  +  EDU184 

aJ 

-  102  (569) 

-  584  (255) 

-  0.000366  (0.00112) 

-  0.00579  (0.00292) 

t,  -  73.8  (43.4) 
t2  -  743  (321) 
w,  -  0.665  (0.172) 

A  -  0.0160  (0.0229) 

EDU1180S  + 

EDU1180R  -t-  EDU184M 

//// 

//// 

EDU1180S  + 

EDU184  •f'  ^UU88SM 

-  24.8  (498) 

-  505  (225) 

-  0.000200  (0.000425) 

-  0.00568  (0.00261) 

t,  -  11.9  (94.6) 
t,  -  462  (170) 
w,  -  0.684  (2.79) 

A  -  0.00616  (0.00479) 

EDU1180S  +  EDU1180R 
-I-  E0U184  -It  EDU885M 


//// 


//// 


APPEHDIX  9 

Optlnized  Parameter  Values  for  the  4-Parameter  Models: 
Any  Category  of  Dive 


EDU885A  + 

EDU885AR  +  EDU885M 


EDU1180S 

+  EDU184  +  EDU885A 
+  EDU885AR  +  EDU885M 
+  DC4W  +  NMR8697 
+  DRREPWET  +  DC8A0W 


model  2 
(2  Tissue 

Mono -exponential) ; 

Pthr  -  0 


T,  -  29.0  (45.2) 

-  623  (179) 

Ai  -  0.00116  (0.000636) 
Aj  -  0.00642  (0.00337) 

-  33.4  (11.4) 

-  715  (108) 

Ai  -  0.00145  (0.000252) 
Aj  -  0.00781  (0.00375) 


model  3 
(1  Tissue 
Bi -exponential) ; 
Pthr  0 


-  14.7  (39.8) 

T2  -  355  (83.1) 

Wi  -  0.839  (0.480) 

A  -  0.00709  (0.00504) 


//// 
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APPENDIX  10 


Log  Likelihoods  and  Optimized  Parameter  Values  for  the  Models 
with  6  or  More  Parameters:  Any  Category  of  Dive 


-ln(likellhood) 


model  4 

model  4 

model  5 

data  set 

null 

o 

1 

k02  0 

k02  -  0 

EDU1180S 

+  EDU184  +  EDU885A 
-»•  EDU885AR  *  EDU885M 
DC4V  -f  NMR8697 
-r  DRREFVET  DC8A0U 

parameter  values: 

394.87 

373.91 

373.41 

376.93 

model  4  ( 3 •mono* exponential ) .  ko2  ■*  0 

T,  -  0.442  (0.174) 

Tj  -  129  (68.9) 
t,  -  767  (162) 

A,  -  0.176  (0.206) 

Aj  -  0.00117  (0.000306) 

-  0.00735  (0.00522) 

model  4  (3*mono«exponential) ,  ko2  *  0 


t,  -  0.284  (0.189) 
tj  -  106  (91.7) 

Tj  -  1180  (599) 

A,  -  0.241  (0.342) 

Aj  -  0.000756  (0.000353) 

A,  -  0.00473  (0.00421) 
ko*  -  0.340  (0.308) 

model  5  (2  bl -exponential) ,  ko2  ~  0 


- 

»A1  - 

A*  “ 


0.125  (6710) 

33.2  (14.0) 

0.256  (10100) 
0.00146  (0.00120) 


Ui 

^■3 

*B1 


^1 


517  (1960) 

2570  (65200) 
0.913  (3.54) 
0.0110  (0.0113) 
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TABLE  1 
Data  Sijomary 


0.7  •  0.7  atB  oxygen,  the  balance  nitrogen 


Mean  valixes 

of  quantities,  averaged  over 

a  data  set,  are 

Indicated  by 

parentheses 

data  set 

DCS  cases 

dives 

gas 

depth 

(fsw) 

bottom 

time 

(min) 

total 

decompression 

time 

(min) 

DCAtf 

18+4(0.5)1/244 

air 

50  -  265 
(154) 

2.9  -  100 
(24) 

3.3  -  99 
(28) 

EDU1180R 

2/128 

0.7 

75  -  151 
(123) 

162  -  270 
(233) 

31  -  176 
(102) 

EOU1180S 

10/120 

0.7 

75  -  150 
(125) 

38  •  126 
(73) 

46  -  176 
(93) 

EDU184 

11/234 

0.7 

40  -  150 
(89) 

20  -  212 
(42) 

2  -  187 
(14) 

EDU88SA 

30/483 

air 

,  50  -  190 
(112) 

14  .  244 
(78) 

1.7  -  290 
(102) 

EDU88SAR 

11/182 

air 

80  -  150 
(102) 

17  -  66 
(N/A) 

2  -  246 
(N/A) 

EDU885M 

4/81 

0.7 

100  -  150 
(133) 

33  -  66 
(51) 

35  .  222 
(76) 

NMR8697 

(11+18(0.5)1/477 

variable 
10-40%  Oj 
the  rest  1 

25  -  130 
.  (69) 

30  -  240 
(112) 

0.6  -  2.5 
(1.5) 

DRREPWET 

3/12 

air 

59  -  177 
(128) 

20  -  40 
(28) 

4.6  -  90 
(50) 

0C8A0V 

(2+0.5)/45  air 

+  O2  dacoopresslon 

90  •  180 
(132) 

2.3  •  60 
(42) 

27  -  106 
(46) 

total 


(92+23(0. 5) J/2006 


TABLE  2 


Log  Likelihoods  for  the  2-,  3*.  and  4*Paraaeter  Models; 
Data  Sets  Not  Cooblned  Across  Categories;  kQ2  ~  0 


-In  (likelihood) 


DCS  . 

cases/dives  null  nodel  1  nodel  2  model  3 


no  Pu  , 


single  air  dives: 

EDU88SA 

+  DC4W  (38+4(0.5)1/727 

single  0.7  dives: 

EDU1180S  *  EOU885M  14/201 

single  dives,  other  gas  mixes: 
NMR8697  (11+18(0.5)1/477 

DC8A0U  (2+1(0.5)1/45 

repeat  air  dives: 

EDU885AK  +  DRREPWET  14/194 

repeat  0.7  dives: 

EDU184  11/234 


154.88 

153.86 

//// 

145.25 

147.70 

50.800 

44.851 

//// 

44.821 

//// 

83.010 

9.655 

81.970 

81.939 

81.844 

81.931 

50.285 

45.719 

44.240 

//// 

44.115 

44.369 

38.455 

37.420 

37.933 

37.329 

52 


TABLE  3 


Log  Likelihoods  for  the  2-,  3-,  and  4-Paraaeter  Models; 
Coobinations  of  Data  Sets  From  Different  Categories 


-In  (likelihood) 


DCS  . 7 . 

cases/dlves  null  model  1  model  2  model  3 


no  w/Pthp 


single,  air  +0.7: 

EDU1180S  +  EOU885A 
+  EDU885M  +  DC4W 

[52+4(0.5)1/928  205.98  204.17  ////  191.21  195.46 

single,  any  gas  mix: 

EOU1180S  +  E0U885A 
+  EDU885M  +  DC4U 
+  DC8A0W  +  NMR8697 

(65+23(0.5)1/1450  299.51  297.23  ////  285.39  //// 


repeat,  air  +0.7: 
EDU184  +  EDU885AR 


+  DRREPWET 

25/428 

95.261 

86.869 

83.544 

85.442 

83.523 

air,  single  +  repeat: 

, 

EDU885A  +  E0U885AR 
+  DC4W  +  DRREPWET 

(52+4(0.5)1/921 

205.55 

202.28 

//// 

194.95 

197.98 

0.7,  single  +  repeat: 

EDU1180S  + 

EDU184  +  EDU885M 

25/435 

101.01 

83.583 

//// 

83.216 

83.395 

any  type  of  dive: 

EDU1180S 

+  EDU184  +  EDU885A 
+  E0U885AR  +  E0U885M 
+  DC4W  +  NMR8697 
■f  DRREFVET 

+  DC8A0W  (90+23(0.5)1/1878  394.87  388.79  ////  377.09  //// 


■ 

TABU  4 

1  Cptlmized  Pmraaeter  Values  for  the  2*  and 

3-Paraaeter  Models: 

Data 

Sets  Not  Cooblned  Across  Categories  I 

■odel  1 

Bodel  1 

(1  Tissue 

(1  Tissue 

Mono-exponential) ; 

Mono  -  exponent  Ip.1  ) 

Pthr  -  0 

Pth«  0 

single  air  dives: 

EDU885A  +  DC4U 

t  -  92.4  (16.8) 

//// 

A  -  0.00222  (0.00035) 

single  0.7  dives: 

EDU1180S  -t-  EOU885M 

t  -  483  (88.6) 

//// 

A  -  0.00586  (0.00260) 

single  dives,  other  gas 

Blxes : 

NMR8697 

t  -  77.9  (34.2) 

t  -  72.9  (28.7) 

A  -  0.00180  (0.000461) 

A  -  0.00269  (0.00414) 

Pth,  -  3.03  (11.2)  1 

DC8A0V 

1 

repeat  air  dives : 

EDU88SAR  •••  DRREPUET 

T  -  297  (83.1) 

T  -  177  (68.7) 

' 

A  -  0.00510  (0.00280) 

A  -  0.0491  (0.0703) 

Pth,  -  6.92  (4.35) 

1  repeat  0.7  dives: 

EDU184 

t  -  413  (90.0) 

t  -  114  (24.9) 

, 

A  -  0.00364  (0.0O194) 

A  -  0.0102  (0.0118) 

Pth,  -  5.34  (3.37)  i 

1 

54 

TABLE  5 


Optlnlzed  Parameter  Values  for  the  2*  and  S-Parameter  Models: 
Combinations  of  Dives  From  Different  Categories 


model  1 
(1  Tissue 

Mono -exponential) ; 

Pthr  -  0 


model  1 
(1  Tissue 
Mono - exponential ) 
Pthr  -‘  0 


single,  air  -1-0.7: 

E0U1180S  +  EDU88SA 
•f  EDU88SM  +  DC4V 


T  -  101  (15.1) 

A  -  0.00223  (0.000307) 


single,  any  gas  mix: 

EDU1180S  -t-  EDU885A 
+  EDU885M  -^  DCAW 
-t-  DC8A0V  -t-  NMR8697 


t  -  99.6  (13.9) 

A  -  0.00214  (0.000255) 


repeat,  air  +0.7: 

EDU134  +  EDU885AR 
+  DRREPWET 


air,  single  +  repeat: 

EDU885A  +  E0U885AR 
+  DC4V  +  DRREPVET 

0.7,  single  +  repeat: 


EDU1180S  + 

E0U184  E0U885M 


a  of  dive: 


EDU1180S 

+  E0Ul84  +  EDU885A 
+  EDU(t85AR  +  EDU835M 
+  DC4W  +  NMR8697 
+  DRREPWET  +  DC8A0W 


T  -  288  (66.2) 

A  -  0.00307  (0.00115) 


t  -  103  (15.8) 

A  -  0.00210  (0.000297) 


T  -  454  (52.9) 

A  -  0.00486  (0.00153) 


t  -  112  (12.7) 

A  -  0.00194  (0.000202) 


t  -  160  (27.9) 

A  -  0.0218  (0.0165) 
fttu  -  6-85  (1.90) 
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TABLB  6 

I 

Optlalzed  Parameter  Values  for  the  4-Farameter  Models: 
Data  Sets  not  Combined  Across  Categories 


miodel  2 
(1  Tissue 

Mono-exponential) ; 

Pthr-  0 


model  3 
(1  Tissue 
Mono - exponent lal ) 
Pthr  0 


single  air  dives: 

EDU885A  -t-  DC4U 

*" 

27.3  (14.6) 

T,  -  14.5  (7.6) 

749  (195) 

T,  -  290  (58) 

Ai  - 

0.00183  (0.00060) 

w,  -  0.894  (  ) 

, 

aJ- 

0.00970  (0.00083) 

A  -  0.00921  (0.00026) 

single  0.7  dives: 

EOU1180S  -I-  EOU885M 

1  1  1  1 

2.51  (310) 

493  (250) 

0.00555  (1.07) 

0.00587  (0.00321) 

//// 

single  dives ,  other 

gas  mixes 

* 

HMR8697 

0.263  (0.955) 

T,  -  1.36  (31.2) 

•»a  “ 

86.7  (54.1) 

T,  -  74.7  (60.7) 

Ai  " 

0.205  (2.57) 

w,  -  0.963  (0.718) 

A2  - 

0.00166  (0.00120) 

A  -  0.00390  (0.0119) 

DC8AOV 

repeat  air  dives: 

EDU88SAR  *  DRREFVET 

//// 

t,  -  141  (92.5) 

T,  -  2900  (7910) 
w,  -  0.925  (0.233) 

A  -  0.160  (0.264) 

repeat  0.7  dives; 

E0U184 

110  (634) 

-  69.3  (54.3) 

*a  “ 

609  (282) 

T,  -  854  (334) 

0.000440  (0.00146) 

Wi  -  0.733  (0.112) 

Aa  “ 

0.00642  (0.00515) 

A  -  0.0237  (0.0374) 

TABLE  7 

Optimized  Parameter  Values  for  the  4- 

Parameter  Models: 

Combinations  of  Data  from  Different  Categories 

model  2 

model  3 

(1  Tissue 

(1  Tlssvie 

Mono -exponential) ; 

Mono - exponent lal ) 

Pth*  -  0 

Pth*  **  0 

single,  air  +0.7: 

EOU1180S  +  EDU885A 

+  EDU88SM  +  DC4V 

■*1  “ 

27.1  (12.6) 

T,  -  13.0  (5.35) 

711  (147) 

Ta  -  312  (39.5) 

Ai  - 

0.00183  (0.000583) 

w,  -  0.900  (0.0461) 

Aj  - 

0.00872  (0.00539) 

A  -  0.0102  (0.0024r) 

\  single,  any  gas  mix: 

EDU1180S  +  EDU885A 

+  EDU885M  +  DC4W 

■*1  "* 

25.6  (10.3) 

//// 

+  DC8A0W  +  NMR8697 

729  (115) 

At  - 

0.00204  (0.000534) 

At  - 

0.00846  (0.00462) 

repeat,  air  +0.7: 

EDU184  +  EDU885AR 

0.304  (0.109) 

T,  -  0.755  (0.813) 

+  DRREPWET 

‘*2  “ 

322  (90.4) 

t;  -  183  (40.1) 

Ai  - 

2.58  (10.8) 

wv  -  0.995  (0.00570) 

Aa  ■■ 

0.00331  (0.00140) 

A  -  0.0351  (0.0332) 

air,  single  +  repeat: 

‘ 

EOU885A  +  EOU885AR 

39.6  (20.2) 

Ti  -  13.2  (7.28) 

+  DC4U  +  DRREPWET 

^2  “ 

732  (212) 

T,  -  249  (63.6) 

Ai  “ 

0.00148  (0.00314) 

w,  -  0.864  (0.0749) 

Aa  " 

0.00807  (0.00677) 

A  -  0.00624  (0.00191) 

0.7,  single  +  repeat: 

EDU1180S  + 

’ 

E0U184  +  E0U88SM 

.;4.8  (498) 

T.  -  11.9  (94.6) 

T, 

505  (225) 

Tj  -  462  (170) 

A  - 

0.000200  (0.000425) 

w,  -  0.684  (2.79) 

Aa  " 

0.00568  (0.00261) 

A  -  0.00616  (0.00479) 

any  type  of  dive: 

EOU1180S 

+  EDU184  +  E0U885A 

■*1  “ 

33.4  (11.4) 

//// 

+  EDU885AR  +  E0U885M 

»2  ■ 

715  (108) 

+  DC4W  +  NMR8697 

Ax  “ 

0.00145  (0.000252) 

+  DRREPWET  +  DC8A0W 

Aa  “ 

0.00781  (0.00375) 

57 

•  '  /  .  ,  - 

TABLE  8 

Predictions  of  P(DCS)  In  Data  Set  EDU1180R  Using  Parameter 
Values  Fitted  to  Other  Data  Sets;  ko2  ~  0 


data  set  to 
which  parameters 
are  fitted 

model 

-LL 

average  P(DCS) 

predicted 

actual 

null 

10.302 

1.56% 

all  single  0.7 

1 

Iw/Pthr 

2 

3 

39.068 

//// 

24.8% 

tl(i' 

//// 

all  repeat  0.7 

1 

Iw/Pth, 

2 

3 

28.409 

32.125 

37.087 

32.029 

17.9 

20.4 

23.5 

20.3 

all  but  EOU1180R 

1 

Iw/Ptju, 

2 

3 

16.199 

MU 

//// 

8.89 

//// 

TABLE  9a 


Predictions  of  P(DCS)  Using  Parameter  Values  Fitted  to  Ither 

Data  Sets;  ko2  -  0 


data  set  to 
which  parameters 
are  fitted 

modal 

data  set  for 
which  predictions 
are  made 

-LL 

averaj'e  P(DCS) 

predicted  actual 

null 

all  repeat  air 

50.285 

7.22% 

all  single  air 

1 

Iw/Pthr 

2 

3 

■  , 

48.788 

52.318 

9.27  (6.37-12.2)% 

(dl  (5.74-12.9) 
11.3  (6.89^15.8) 

null 

all  repeat  0.7 
except  EDU1180R 

44.369 

4.70% 

all  single  0.7 

1 

Iw/Pthr 

3 

38.936 

3^7 

//// 

5.28  (2.43-8.14)% 

(0-100) 

//// 

null 

all  single  0.7 

50 . 800 

6.97% 

all  single  air 

1 

Iw/Ptte 

2 

3 

50.496 

48.703 

6.82  (4.59-9.06) 

(4.05-8.32) 
5.58  (3.27-7.88) 

null 

all  repeat  0.7 
except  EDU1180R 

44.369 

4.70% 

all  repeat  air 

1 

Iw/Pthr 

2 

3 

44.299 

43.124 

//// 

45.254 

10.6  (4.96-16.3)% 
10.3  (0.79-19.9) 

lo{8'^(2.23-19.5) 

■ 

null 

all  repeats 
except  EOU1180R 

86.080 

5.84% 

all  singles 

1 

92.634 

8.64  (6.66-10.6)% 

Iw/Pthr 

3 

//// 

(6.76-12.0) 

//// 

null 

all  0.7  dives 
except  EOU1180R 

101.01 

5.75% 

all  air  dives 

Iw/Ptta 

3 

93.075 

94.971 

7.01  (5.58-9.58)% 

(5.14-8.74) 
8.35  (6.16-10.5) 

TABLE  9b 


Predictions  of  P<DCS)  Using  Parsmeter  Values  Fitted  to  Other 
Data  Sets;  ^  -  0  --  continued 


Note:  The  ln( likelihood)  is  undefined  when  the  model  predicts  zero, 
risk  for  a  dive  in  which  a  case  of  DCS  was  recorded,  because  in  this 
sitxiation  the  likelihood  of  the  observed  family  of  dive  outcomes, 
given  that  P(DCS)  is  governed  by  the  model  in  question,  is  zero. 


data  set  to 
which  parameters 
are  fitted 

model 

data  set  for 
which  predictions 
are  made 

-LL 

average  F(DCS) 

predicted  actual 

null 

all  single  air 

154.88 

5.50% 

all  repeat  air 

1 

Iw/Pthr 

3 

undefined 

undefined 

//// 

undefined 

7.23  (2.81-11.7)% 
11.1  (2.65-19.3) 
//// 

9.90  (?.50-17.2) 

null 

all  single  0.7 

50.800 

6.97% 

all  repeat  0.7 
except  EDU1180R 

Jw/Ptht 

2 

3 

45.504 

50.046 

45.949 

undefined 

5.77  (2.40-9.13)% 
8.18  (0.32-16.1) 
5.94  (1.73-10.2) 
4.31  (0-10.4) 

null 

all  single  air 

154.88 

5.50% 

all  single  0.7 

3 

ijndefined 

//// 

151.42 

//// 

4.48  (2.29-6.69)% 

5^^2^'^(0-38.8) 

//// 

null 

all  repeat  air 

50.285 

7.22% 

all  repeat  0.7 
except  EDU1180R 

Iw/Pthr 

3 

' 

52.829 

48.061 

52.532 

51.207 

2.73  (1.15-4.29)% 
3.37  (0.66-6.13) 
3.82  (0.46-7.18) 
3.32  (1.18-5.47) 

null 

all  singles 

299.51 

5.28% 

all  repeats 
except  EDU1180R 

Iw/Pth* 

3 

undefined 

undefined 

319.88 

undefined 

4.61  (2.72-6.50)% 
7.74  (3.15-12.3) 
7.49  (0-22.7) 

7.36  (2.96-11.8) 

null 

all  air  dives 

205.55 

5.86% 

all  0.7  dives 
except  EDU1180R 

Iw/Pthr 

3 

undefined 

Mix 

undefined 

3.89  (2.42-5.38)% 

4^{^^(0.61-7.74) 
3.92  (2.38-5.48) 

TABLE  10 


Log  Likelihoods  and  Optimized  Parameter  Values  for  the  Models 
with  6  or  More  Parameters:  Any  Category  of  Dive 


data  set 
EDU1180S 

•f  EDU184  +  EDU885A 
+  EDU885AR  +  EDU88SM 
+  DC4W  +  NMR3697  . 

+  DRREPWET  +  DC8A0V 


-ln(likellhood) 


null 

394.87 


model  4 

ko2  -  0 


model  4 

ko2  "  0 


model  5 

ko2  -  0 


373.91 


373.41 


376.93 


parameter  values: 

model  4  (3-mono-exponentlal) ,  ko2  -  0 

T,  -  0.442  (0.174) 

-ri  -  129  (68.9) 
ta  -  767  (162) 

-  0.176  (0.206) 

a:  -  0.00117  (0.000306) 

-  0.00735  (0.00522) 

model  4  ( 3 -mono 'exponent lal) ,  kga  *  0 
Ti  -  0,284  (0.189) 


^2  •• 
T,  •« 


106  (91.7) 
1180  (599) 


a:  -  0.241  (0.342) 

Aa  -  0.000756  (0.000353) 
a1  -  0.00473  (0.00421) 
l4  -  0.339  (0.308) 

model  5  (2  hi - exponent ial ) ,  koj  “  0 

tA,  -  0.125  (6710) 
tA;  -  33.2  (14i0) 
wAi  -  0.256  (10100) 

Aa  -  0.00146  (0.00120) 

T-i  -  517  (1960) 
tS:  -  2570  (65200) 

-  0.213  (3.54) 

Ab  -  0.0110  (0.0113) 
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